The asymptotic iteration method for the angular spheroidal eigenvalues with 

arbitrary complex size parameter c 



T. Barakat* K. Abodayeh, B. Abdallah, and 0. M. Al-Dossary 
Physics Depntartment, King Saud University 
Riyadh, Saudi Arabia 



Abstract 

The asymptotic iteration method is applied, to calculate the angular spheroidal 
eigenvalues A™(c) with arbitrary complex size parameter c. It is shown that, 
the obtained numerical results of A™(c) are all in excellent agreement with 
the available published data over the full range of parameter values £, m, and 
c. Some representative values of A™(c) for large real c are also given. 
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1 Introduction 



The solution of the spheroidal wave equation is a very old subject, but it still an 
important theme in the existing literature. The importance of this equation arises 
in many areas of physics. For instance, it plays an important role in the study of 
light scattering in optics [1-3], nuclear modeling [4], signal processing and commu- 
nication theory [5], electromagnetic modeling [6], and in finding the electromagnetic 
induction (EMI) response of canonical objects at magnetoquasistatic frequencies [7]. 
Applications utilizing complex c include for example, light scattering from spheroidal 
particles, and spheroidal antennas enveloped in a plasma medium. 

Attempts to find rapid, and accurate eigenvalues A™(c) of the angular spheroidal 
wave equation for large size parameter c 2 (assumed real) have been ongoing. Flam- 
mer summarizes the work (up to 1957 [8]), and documents the asymptotic expansions 
for A™(c) of the angular spheroidal wave equation. Since that time, serious attempts 
for this case were made by many authors. Slepian [9], and Streifer [10] derived uni- 
form asymptotic expansions for the spheroidal functions and their eigenvalues, which 
were further developed by des Cloiseaux, and Mehta [11], and Dunster [12]. Other 
asymptotic results based on WKB methods have been obtained by Sink, and Eu 
[13]. Recently, asymptotic expansions of A™(c) for large /, and c have been proposed 
by Guimaraes [14], de Moraes and Guimaraes [15], and the work of Do-Nhat [16, 
17] summarizes and provides more details of Flammer's expansions for A™(c). 

Nevertheless, because of the complexity of the angular spheroidal wave equation 
with arbitrary complex size parameter c = c r + Cji, where c r = Re{c}, and q = 
Im{c}, evaluation of A™(c) in this regime has been much less studied. 

Very recently, Barrowes et al. [18] compute the asymptotic expansions of A™(c) 
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with arbitrary complex size parameter c in the asymptotic regime of large | c | 
with /, and m fixed. On the other hand, few packages have been developed for the 
computation of the angular spheroidal eigenvalues A™(c) with arbitrary complex size 
parameter c = c r + Cii. Thompson [19], Li et al. [20, 21], and Falloon et al. [22] are 
of the most recent ones. 

Those attempts to obtain the eigenvalues A™(c) with arbitrary complex size 
parameter c rely heavily on power series expansions, and complicated 

recurrence relations. Accurate results in those works are obtainable at the expense of 
extensive mathematical, and numerical manipulations, thus obscuring the physical 
analysis of the corresponding system. 

The present work applies the asymptotic iteration method (AIM) [23, 24], for the 
computation of the angular spheroidal eigenvalues A™(c) with arbitrary complex size 
parameter c = c r + c^. This method was applied by Barakat et al. [25] to compute 
the angular spheroidal eigenvalues A™( c ) with real c 2 , and for the eigenenergies of 
the anharmonic oscillator potential [26]. The implementation of this method was 
straightforward, and the results were sufficiently accurate for practical purposes. 
Most importantly, the numerical computation of the angular spheroidal eigenvalues 
using this method was quite simple, fast, and the eigenvalues were satisfying a simple 
ordering relation. Therefore, one can unambiguously select the correct starting 
eigenvalue. 

Furthermore, AIM was quite flexible in the sense that, it is applicable to any 
parameter value involved like £, m, and c. It also handles A™(c) with large £, and 
c which poses many numerical instabilities to some of the previously mentioned 
methods. Therefore, the main motivation of the present work is to overcome the 
shortcomings of those approaches, and to formulate an elegant algebraic approach 
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to yield a fairly simple analytic formula which will give rapidly the eigenvalues with 
high accuracy. 

In this spirit, this paper is organized as follows. In Sec. 2 the asymptotic 
iteration method for the angular spheroidal wave equation is outlined. The analytical 
expressions for asymptotic iteration method are cast in such a way that allows the 
reader to use them without proceeding into their derivation. In Sec. 3 we present 
our numerical results compared with other works, and then we conclude and remark 
therein. 

2 Formalism of the asymptotic iteration method for the angular spheroidal 
wave equation 

The angular spheroidal wave equation, with which we shall be concerned, is 

Se, m (c; rj) = 0. (1) 

The parameter c, which is related to the ellipticity of the spheroidal coordinate 
surfaces, is allowed to be a complex variable in the present work. Consequently, 
there is no need to make the usual distinction between the prolate, and oblate forms 
of the spheroidal wave equation, and the prolate form equation (1) is adopted for 
definiteness. The other parameters A™(c), and m are separation constants. 

The second arises as a wave number for the polar angle of spheroidal coordi- 
nates and, as usual, is required to be a nonnegative integer; t > m is an integer 
enumerating the eigenvalues, and functions. 

The spheroidal wave functions Se jTn (c; rf) are defined to be the solutions of equa- 
tion (1) that are finite at the two end points rj — ±1 of the range of the independent 
variable. These finiteness can be satisfied only for certain eigenvalues A™(c), which 
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depend on the values of c once a specific value of m has been chosen. 

The simplest case is that of c = 0, for which the function S^ m (c; 77) reduces to 
the associated Legendre function, and (A™(c)) 2 = £(£ + 1) is its eigenvalues. 

Here, the integer i labels successive eigenvalues for fixed m. When I = m we 
have the lowest eigenvalue, and the corresponding eigenfunction has no nodes in 
the interval — 1 < 77 < 1. When I = m + 1 we have the next eigenvalue, and the 
eigenfunction has one node inside (-1, 1); and so on. A similar situation holds for 
the general case 0. 

In order to apply the AIM for the general case c 7^ 0, we have to investigate 
the behavior of the solution near the singular points 77 = d=l. Substituting a power 
series expansion of the form 

00 

s e , m (c;v) = (i-v 2 Tj2 a ^ 1 - r i 2 )\ (2) 

into equation (1), we find that the regular solution has a = m/2. Without loss 
of generality, we can take m > since m — > — m is a symmetry of the equation. 
Therefore, we get an equation that is more tractable to the method if we factor out 
this behavior. Accordingly, we set 

s e , m (c;v) = (i-v 2 ) m/2 yUc;v), (3) 

then the new function ^, m (c; r/) will satisfy a second-order homogenous linear dif- 
ferential equation of the form 

(1 — ^7 ) ^2 2 ( m + l W ^ + {e-cr] )ye, m {c; r]) = 0, (4) 

where 

e = (Xl\c)f-m{m + l). (5) 
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Both equations (1), and (4) are invariant under the replacement rj — > —r\. Thus 
the functions Se jm (c; rf), and ye, m (c; rf) must also be invariant, except possibly for an 
overall scale factor. 

The systematic procedure of the AIM begins by rewriting equation (4) in the 
following form 

y'e, m (^ v) = Mv)y't,m(c; v) + so(ri)ye,m( c ; v), ( 6 ) 

where 

2(m + 1>7 e - c 2 if 

X °^) = (1,^2) » and s o(v) = ~ ( i_ r] 2y ( 7 ) 

and following the technique of AIM [23, 25], that will lead to a general solution of 
equation (6): 

ye, m (c; r]) = exp ^- J (3(rj)drj^j C 2 + C x j exp ^ J {X (r)") + 2(3(r)")}dr)"^J drj' 
If for some n > 0, 

0(v) = = Sn - M (9) 

with 

^n(v) = Ki-iiv) + s n -i(v) + A (?7)A n _i(?7), and s„(r]) = s' n ^(r]) + s (?7)A„_i(?7).(10) 

For sufficiently large n, we can now introduce the termination condition of the 
method, which in turn, yields the angular spheroidal eigenvalues A™(c) 

S n (v) = s n (v)^n-i(r]) ~ s n -i(r])K{r]) = 0- (11) 



(8) 
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3 Numerical results for the angular spheroidal eigenvalues A™(c) 

Within the framework of the AIM mentioned in the above section, the angular 
spheroidal eigenvalues A™(c) are calculated by means of equation (11). To obtain 
the eigenvalues A™(c), the iterations should be terminated by imposing a condition 
dn{v) = as an approximation to equation (11). On the other hand, for each itera- 
tion, the expression 5 n (r)) = s n (r))\ n -i(r)) — s n -i(r))\ n (r)) depends on two variables: 
A™(c), and rj. The calculated eigenvalues A™(c) by means of this condition should, 
however, be independent of the choice of rj. Nevertheless, the choice of rj is observed 
to be critical only to the speed of the convergence to the eigenvalues, as well as for 
the stability of the process. In this work it is observed that, the best starting value 
for f] is the value at which the effective potential of equation (1) takes its minimum 
value. For this purpose, it is necessary to perform the variable change 77 — > tanh(x), 
mapping the finite interval (-1,1) into the infinite one (—00, 00), and then equation 
(1) can be rewritten as 

d 2 

- fa2 S *,™ + v eff( x ) s e,m = -m 2 S^ mi (12) 

where the effective potential V e ff(x) is 

V eff (x) = - [{\T{c)) 2 - c 2 ] sech 2 {x) - c 2 sech 4 (x). (13) 

V e ff(x) is an even function, its minimum value occurs when x — 0, which in turn 
implies that i] = 0. Therefore, at the end of the iterations we put rj = 0. 

To test the rate of convergence of AIM numerically, we calculate the angular 
spheroidal eigenvalue Aq(10) shown in figure 1. We simply tried a set of iteration 

numbers n — 5, 10, , and the convergence of AIM seems to take place smoothly 

when n > 45 iterations. 
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Proceeding in the same way, the results of the AIM for A™(c) with different values 
of £, m, and c are reported in tables 1, 2, 3, 4, and 5. The results are shown in such 
a way that one can easily judge the accuracy of the method. The angular spheroidal 
eigenvalues A™(c) were calculated by means of 45 iterations only. Our calculated 
eigenvalues A™(c) are all in excellent agreement with the available published data 
[8, 18, 20, 22] over the full range of parameter values, £, m, and c. 

A second, more stringent, test of the method is shown in Table 3, where we 
consider large real c, in this case, in order to reproduce more accurate results, the 
angular spheroidal eigenvalues A™(c) were calculated by means of 100 iterations. 
Again the agreement is quite excellent with the available published data. 

It is worthwhile to emphasize that, the AIM is very easy to implement for cal- 
culating the angular spheroidal eigenvalues A™(c), without having to worry about 
the ranges, and forms of c. Moreover, for purely real c, and purely imaginary c, the 
obtained eigenvalues are satisfying a very simple ordering relation. But, for complex 
arbitray c values, the eigenvalues are no longer real, since the spheroidal equation is 
not self-adjoint. In this case the ordering is determined by the absolute magnitude of 
(A™(c)) 2 . Hence, one can unambiguously select the correct starting eigenvalue. This 
represents a significant advantage over the tridiagonal matrix method [22] in which, 
the eigenvalues are not ordered, and hence to choose the correct matrix eigenvalue, 
one must use an iterative process to move towards the starting value. 

As a concluding remark, we would like to point out that, the accuracy of the 
results could be increased if the number of iterations are increased. 

This work is supported by King Saud University, College of Science - Research 
center projects No. (PHYS/2005/26). 
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Table 1: Comparison of selected values of eigenvalues (A™(c)) 2 computed by Flam- 
mer [8], Le-Wei Li et al. [20], and by means of the present work. 













c 2 




Flammer [8] 


Le-Wei Li et al. [20] 


Present work 


-1.00 


(4,11) 


131.560 


131.560 


131.560 


0.10 


(2,2) 


6.01427 


6.01427 


6.01427 


1.00 


(1,1) 


2.19555 


2.19555 


2.19555 




(2,2) 


6.14095 


6.14095 


6.14095 




(2,5) 


30.4361 


30.4362 


30.4362 


4.00 


(1,1) 


2.73411 


2.73411 


2.73411 




(2,2) 


6.54250 


6.54250 


6.54250 


16.00 


(1,1) 


4.39959 


4.39959 


4.39959 




(2,5) 


36.9963 


36.9963 


36.9963 
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Table 2: Comparison of selected values of eigenvalues (A™(c)) 2 computed by Falloon 
et al. [22], and by means of the present work. 





c 


(m,£) 


Falloon et al. [22 


Present work 


10 


(0,0) 


9.228304 


9.228304 




(0,1) 


28.13346 


28.13346 




(0,2) 


45.86895 


45.86895 




(1,1) 


10.28777 


10.28777 




(1,2) 


29.33892 


29.33892 




(1,3) 


47.30152 


47.30152 




(2,2) 


13.46308 


13.46308 




(2,3) 


32.93818 


32.93818 




(2,4) 


51.52485 


51.52485 


lOi 


(0,0) 


-81.02794 


-81.02794 




(0,1) 


-81.02794 


-81.02794 




(0,2) 


-45.48968 


-45.48968 




(1,1) 


-62.11935 


-62.11935 




(1,2) 


-62.11915 


-62.11915 




(1,3) 


-29.18576 


-29.18576 




(2,2) 


-43.29025 


-43.29025 




(2,3) 


-43.28716 


-43.28716 




(2,4) 


-13.50811 


-13.50811 
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Table 3: Comparison of selected values of eigenvalues (A™(c)) 2 computed by J. 
Liu [27], Falloon et al. [22], and by means of the present work. 



c 


(m,£) 


J. W. Liu [26 


Falloon et al. [22 


Present work 


50 


(0,0) 


49.24615 




49.24615 




(0,1) 


148.2306 




148.2306 




(0,3) 


343.1109 




343.1109 




(0,4) 


438.9725 




438.9725 




(1,1) 


50.25646 




50.25646 




(1,2) 


149.2622 




149.2622 




(1,4) 


344.1894 




344.1894 




(1,5) 


440.0769 




440.0769 


100 


(0,0) 


99.24810 


99.24810 


99.24810 




(0,1) 


298.2405 


298.2405 


298.2405 




(0,2) 




496.2212 


496.2212 




(0,3) 


693.1825 




693.1825 




(0,4) 


889.1162- 




889.1162 




(1,1) 


100.2532 


100.2532 


100.2532 




(1,2) 


299.2558 


299.2558 


299.2558 




(1,3) 




497.2472 


497.2472 




(1,4) 


694.2195 




694.2195 




(1,5) 


890.1645 




890.1645 
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Table 4: Comparison of selected values of eigenvalues (A™(c)) 2 computed by Le-Wei 
Li et al. [20], and by means of the present work. 



me)? 

m,f) Le-Wei Li et al. [20] Present work 



1.824770+2 
2.094267+5 
5.217093+3 
3.563644+2 
1.998555+4 
3.862833+4 
2.136987+5 



.601670i 
.807965i 
.081362i 
.887165i 
.097453i 
.492300i 
.449457i 



0,0) 1.701836+4.219998i 

0,2) 1.993901+8.576325i 

0,2) 23.91023+18.74194i 

0,1) 10.13705+11. 12216i 

1,1) 2.919098+6. 134851i 

1,1) 12.19691+16.24534i 

2,0) 6.098946+7.684379i 



1.701836+4.219998i 
1.993901+8.576325? 
23.91033+18. 74184i 
10. 13705+11. 12218? 
2.919095+6.13485H 
12.19691+16.24534i 
6.098961+7.684333i 
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Table 5: Comparison of selected values of eigenvalues (A™(c)) 2 computed by B. E. 
Barrowes et al. [18], and by means of the present work. 



m,£ 



B. E. Barrowes et al. [18] Present work 



1.824770749208805+ 2. 
3.563644553545243+ 2. 
5.217093042404772+ 3. 
4.067274712533398+ 6. 
2.244329796261236+ 8 
7.606334073445308+ 6 
6.316233767329015+ 9 



601670692890318*- 
887165344336900*' 
081362886557631* 
264358978587767* 
973752190228394* 
906465157219409* 
949229739353585*' 



0,0 
0,1 
0,2 
0,3 
0,4 
0,5 
0,6 



I. 705180+4.220186* 
10.14084+11.12159* 
23.91583+18.74332* 

II. 78093+22.54139* 
2.156125+12.81092* 
48.80665+54.01199*- 
29.84604+55.68999* 



I. 705180+4.220186* 
10.14084+11.12159* 
23.91583+18.74332* 

II. 78093+22.54139* 
2.156124+12.81092* 
48.80665+54.01197*- 
29.84604+55.68998*' 



1.998555442181652+ 4. 
3.862833529248772+ 4. 
2.184204069300826+ 7. 
7.270040170458184+ 5. 
6.119087892218941+ 8. 
4.510843794687041+11 



097453662365392*- 
492300074953849* 
326156812534641* 
010809182556227* 
234638882858787*' 
068777156965684*- 



1,0 

1,1 
1,2 
1,3 
1,4 
1,5 



2.915319+6.133951*' 
12.20110+16.24408* 
3.102506+10.53921*- 
47.41099+39.42618* 
29.88062+46.18128*- 
14.19982+38.58584* 



2.915319+6.133951*' 
12.20110+16.24408* 
3.102506+10.53921* 
47.41098+39.42619*- 
29.88062+46.18128* 
14.19982+38.58584* 



2.136987377094029+ 5. 
4.105156484215650+ 5. 
5.907125751703487+ 6. 
7.634658702877481+ 6. 
6.337223309080594+ 9. 



449457313914277* 
922750658440496* 
283464345814330*- 
576768822063829* 
739915760209533*- 



2,0 

2,1 
2,2 
2,3 
2,4 



6.102540+7.684764*- 
16.13687+20.40623* 
32.08759+34.48753* 
53.56137+49.63852*- 
34.18267+53.57326*- 



6.102540+7.684763*- 
16.13687+20.40623*- 
32.08758+34.48754*- 
53.56136+49.63853*- 
34.18267+53.57326*- 



2.254441944326160+ 6. 
4.312789375877335+ 7. 
6.178462421804212+ 7. 
7.954733051349589+ 8. 



731940814252908*- 
267942971679416* 
686866389128842* 
033437876318695* 



3,0 
3,1 
3,2 



11.27374+9.046369*- 
21.98994+24.06425* 
38.86809+40.69989* 
61.43614+58.62467* 



11.27374+9.046369*- 
21.98994+24.06425* 
38.86808+40.69989* 
61.43613+58.62466* 



2.357663561227191+ 7 
4.496463440238013+ 8. 
6.420513778117898+ 9. 



971913317957412* 
560827522001995* 
029346203504799* 



4,0 
4,1 
4,2 



18.43306+10. 
29.78282+27. 
47.52324+46. 



28552*- 
39486* 
358 17i 



18.43306+10. 
29.78282+27. 
47.52324+46. 



28552*- 
39486* 
35817*- 



2.450444507315414+ 9. 
4.662342487692823+ 9. 



182664291501323*- 
817718575239088*- 



5,0 
5,1 



27.58314+11. 
39.52918+30. 



43648*- 
48922i 



27.58314+11. 
39.52918+30. 



43648*- 
48922i 



2.535162563188484+10.371846133322299i 



6,0 



38.72574+12.51971*- 



38.72574+12.51971*- 
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n 



Figure 1: The rate of convergence of AIM for Ag(10) as a function of the number of 
iterations n. 
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